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We present a general derivation of Hess-Fairbank effect or nonclassical rotational inertial (NCRI) , 
i.e. the refusal to rotate with its container, as well as the quantization of angular momentum, as 
consequences of off-diagonal long-range order (ODLRO) in an interacting Bose system. Afterwards, 
the path integral formulation of superfluid density is rederived without ignoring the centrifugal 
potential. Finally and in particular, for a class of variational wavefunctions used for solid helium, 
treating the constraint of single-valuedness boundary condition carefully, we show that there is no 
ODLRO and, especially, demonstrate explicitly that NCRI cannot be possessed in absence of defects, 
even though there exist zero-point motion and exchange effect. 
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I. INTRODUCTION 

It was first suggested by London that the ability of 
liquid ^He II to flow through narrow capillaries without 
apparent friction is a consequence of Bose-Einstein con- 
densation (BEC)ii The concept of BEC was later general- 
ized by Penrose and Onsager to be applicable to interact- 
ing particles. ^'"^ It was further generalized and systemat- 
ically investigated by Yang, as the notion of off-diagonal 
long-range order (ODLRO):- Now it is known that the 
no-friction behavior in narrow capillaries is only one of 
several phenomena of superfluidityi^ As elaborated by 
Leggett,^ the most basic manifestation of superfluidity is 
the Hess-Fairbank effect,^ which was also called "nonclas- 
sical rotational inertial" (NCRI) by Leggetti^ This refers 
to the refusal of the system to rotate with its container, 
when its angular velocity is sufficiently low. It is the 
counterpart of the Meissner effect of superconductivity. 
Furthermore, the quantization of angular momentum of 
the superfluid in the rotating container is the counterpart 
of the magnetic flux quantization in a superconductor. 

In the case of superconductivity, the demonstration of 
Meissner effect and the magnetic flux quantization, as 
consequences of ODLRO, was made by Yangj4 and by 
Sewell and Nieh et al. in a more recent alternative ap- 
proachiS Bloch discussed the relation between supercon- 
ducting persistent current and ODLRO— In the case of 
superfluidity, Kohn and Sherrington derived the Hess- 
Fairbank effect as a consequence of ODLRO by using 
a sophisticated hierarchy of equations of thermal Green 
functions. For a noninteracting Bose gas and a Gross- 
Pitaevskii system, Leggett made a clear-cut demonstra- 
tion of Hess-Fairbank effect and quantization of angu- 
lar momentum as consequences of BEC^ Earlier, in an 
extremely thorough and insightful discussion,^ Leggett 
pointed out that a sufficient condition of superfluidity is a 
certain topological connectedness property of the many- 
body wavefunction, and that at least for zero tempera- 
ture, ODLRO gives rise to this connectivity and thus su- 
perfluidity, but for a finite temperature, whether ODLRO 
is sufficient for superfluidity in general is not conclusive. 



Moreover, Leggett established, from the point of view 
of connectivity of wavefunction, that BEC and NCRI be- 
havior can in principle also be exhibited by a solid^ 
Recently, Kim and Chan clearly observed NCRI-like be- 
havior in bulk solid ^He in an annulus channcl,-'^'^ shortly 
after an earlier such observation in solid ^He confined in 
porous Vycor glassii^ But a consensus on its origin is yet 
to be reached The earliest predictions on super- 
solidity, i.e. superfluid behavior in a solid, were based 
on BEC of defect statesii^iSS But the concentration of 
zero-point vacancies is less than 10^^ according to the 
experimental resultsi^ Thus an important question is 
whether it is possible for a pure commensurate sample 
of solid *He, i.e. without vacancies or interstitials, be- 
come a supersolid. Negative answers were given recently 
in a path integral Monte Carlo calculation of exchange 
frequencies in bulk hep '^He^ and in a general argument 
about superfluidity densityi^i 

Thus from both the fundamental point of view and 
the perspective of understanding supersolid behavior, it 
appears still interesting to make a general derivation of 
the Hess-Fairbank effect and quantization of angular mo- 
mentum as clear consequences of ODLRO for an inter- 
acting Bose system in a rotating container. In this pa- 
per, we first make such a derivation. Afterwards, for a 
reason explained below, we rederive the superfluid den- 
sity in the path integral formulation which is the very 
basis of the analyses of solid ^He in Refi and Ref,— . 
Finally, we consider the trial wavefunctions ever used in 
variational calculations for solid helium, including the 
Hartree wavefunction, the Hatree-Fock wavefunction and 
the Nosanow-Jastrow wavefunction. It is shown that 
there is no ODLRO or BEC in these wavefunctions. 
Moreover, by examining the dependence of free energy 
on the rotation velocity of the container, we explicitly 
demonstrate that a commensurate solid described by 
such wavefunctions cannot possess NCRI, in absence of 
vacancies or interstitials, even if there exist zero-point 
motion and the exchange effect. 

Note that the non-superfluidity of the Hartree- 
Fock wavefunction made up of localized single parti- 



cle wavepackets has been discussed by Leggett from 
the point of view of disconnectivity of the wavefunction 
long agoi^ Our approach provides an exphcit construc- 
tion of the rotating wavefunction under the constraint of 
the "single-valuedness boundary condition" (SVBC) as 
called by Leggett^, while keeping the energy the same 
as that in the static case. To do this, adjustment on 
the wavefunction needs to be made in the exponentially 
vanishing regions, indeed as argued by Leggett. 

The organization of the paper can be clearly seen in 
the section titles. 



II. HAMILTONIANS AND FREE ENERGIES IN 
THE TWO REFERENCE FRAMES 
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FIG. 1: The cylindrical annular container as often considered 
in literature and also here. The radius R is much larger than 
the thickness d. The rotation is along the axis of the two 
concentric cylinders. 



As usual, consider a Bose system in a container rotat- 
ing with a angular velocity a;. Thermodynamic equilib- 
rium is determined by the minimization of the free energy 
in the co-rotating frame of reference, in which the wall of 
the container is at rest. In this frame, the Hamiltonian 



H 



(pj — mu: X Vj)'^ 1 



2m 



(1) 

where the notations are standard, U is the external po- 
tential, Vjk = V{\rj — r^l) is the particle-particle inter- 
action and is rotationally invariant. For basic mechanics 
and thermodynamics of a rotating body and the applica- 
tion to a Bose system, we refer to the standard textsiS^ 
But we draw attention to the point that for each parti- 
cle, the radius vector r^, the canonical momentum Pj, 
and the angular momentum 1^ = Vj x are respectively 
the same in the laboratory frame and in the co-rotating 
frame. It is for this reason that H can be re-written as 



H — Hiab — <^ ■ ^ Ij 



where 



ff/a6 = 5][p,V2m + C/(r,)] 



is the Hamiltonian in the laboratory frame. This point 
is quite delicate in ODLRO studyiSi 

For simplicity, as usual, consider a thin cylindrical 
annular container, with average radius R and thickness 
d R (Fig. The rotation uj is, of course, along the 
cylindrical axis (z axis). Then the centrifugal potential 
becomes a w-dependent constant (in the sense that it is 
independent of the particle configuration), — iM(tji?)^, 
where M is the total mass of the particles. 

It is probably useful to make a synopsis here on the free 
energies in the two reference frames and their relations 
with the rotational inertial and the superfluid density. 



The free energy in the co-rotating frame can be written 
as 



(2) 



where a is a constant, Ic — MR^ is the classical rotation 
inertial, 

Fo = a+^{Ic-I)LJ^. 
The total angular momentum is L = Ij), hence 



L-, = Ilu = — 



dF 

dijj 



In the laboratory frame, the free energy is 
Fiab^F + uj-'L = F + Ioj^. 

Therefore, 



where 



1 9 1 9 

Flab = Fiabfi + 2^cW = « + 2^"^ ' 



Flabfi = a - -j^ilc - 



Consistently, one also has 



L, ^ luj = 



dFiab 



I = - 



d^F d^Fiab 



For a normal system, / = Ic, thus F^ = Fiabfl = a- If 
Fq or, equivalently, Fiabfl depends on lo, then the system 
is a superfluid, with NCRI. The superfluid fraction is 

PS _^ I _ Id^Fp _ 1 d^Fiabfl 
p Ic Ic duj^ Ic du!^ ' 

where ps and p are the superfluid density and the total 
fluid density, respectively. 

It should be noted that in equilibrium, it is F, not 
Flab, that is related to the partition function Q as F ^ 
-kTlnQ. 
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III. A DERIVATION OF NCRI FROM ODLRO 

Now we make a general derivation that if the system 
possesses ODLRO, then Fq in Eq. depends on oj. We 
use an approach similar to Yang's treatment of supercon- 
ductivity in a magnetic field^ 

Using the cylindrical coordinates (z, r, 9) and consider- 
ing the geometry, the Hamiltonian Q can be simplified 
as 



(3) 



with 



{pej - mLoRf ^ Pij^ ^ 1 



2m 



2m 



where pgj = {\/ R)d / dOj. The radial momentum prj = 
d/drj is neglected because d <^ R. An eigenfunction 
of H satisfies 

Hip a = Ealpa, 

and the periodic boundary condition, or SVBC, 

+ 2^, {9,^,}) = i^o.{6,,{0^^J}) (4) 

due to the cylindrical geometry. 

Because —MuP'R? 12 is a constant for a given u), we 
only need to consider Hq^ whose eigenfunctions are com- 
pletely the same as those of , i.e. 



where = + Muj'^R^/2. 
By a "gauge" transformation 



(5) 



V-a = -(/-a exp( 



imujR Y,j 6] . 



^'^ satisfies 



a' 



where 



\Y.^,k+NU{R). 



2m 2m' 2 
The angular boundary condition becomes 

+ 27r, {9.^,}) = e-'^M9,, {9.^,}). (6) 
Now consider the un-normalized density matrix 



Pdn 



Because H = Hq ~ IcW^/2, Fq in Eq. Q is given by 
Fq = -kTlnQ{Ho), where Q{Ho) = Trpdm- From pdm, 



by tracing over all but one particle, one obtains the (un- 
normalized) one-particle reduced density matrix pi. 

The problem determined by Hq together with SVBC 
is equivalent to the description in terms of H'q together 
with Eq. ||HJ). From the cj-independence of i/p and the 
boundary condition ©, one knows that 



27r|pi| 



9'|pi|6l-27r) =e" 



'-{9'\p,\9). (7) 



We can now apply Yang's method to the current prob- 
lem. Without ODLRO, pi is vanishingly small except in 
the regions around 9 = 9' zL 2mT, where n = 0, 1, • • •. As 
indicated by Eq. {Tj) , the values of pi in two neighboring 
regions only differ by a phase factor ""r" . 

With ODLRO, these regions with nonvanishing pi 
merge into each other, and pi is nonvanishing every- 
where. The above phase relation remains. 

Furthermore, with ODLRO, Eq. ([TJl implies that the 
dependence of (r'|pi|r) on r — r' must vary as lu varies. 
Consequently, Q{Ho) and thus _Fb also vary with lu, as 
Q{Ho) = Tripi and Fq = - kT In Q{Ho). This proves 
that ODLRO gives rise to superfluidity or NCRI. 



IV. QUANTIZATION OF ANGULAR 
MOMENTUM 

We now demonstrate the quantization of angular mo- 
mentum as a consequence of ODLRO, by employing the 
method of Bloch in discussing superconducting persistent 
current,^ and also as a generalization of an argument by 
Leggettfli As said above, the angular momentum and 
momentum are, respectively, the same in the laboratory 
frame and in the co-rotating frame. But for convenience, 
here we use the co-rotating frame. 

Consider the one-particle reduced density matrix with 
r and z coordinates integrated over. 



{9'\pi\9) = J drj dz{r,9',z\pi\r,9,z), 
and its Fourier transformation 



where I and I' represent angular momenta. Its normal- 
ization is 



{9\p,\9)d9^J2(l\p,\l)^N. 



Conservation of angular momentum in the z direction 
implies that {l'\pi\l) = for I' ^ I. 

The total angular momentum, along the z direction, 
for the system under consideration can be given as 



L, = J2l{l\p^\l). 



(9) 
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In the Hamiltonian in Eq. (PJ, pgj can be substituted 
as l^j/R, where l^j is the z-component angular momen- 
tum operator of the j-th particle. Thus H depends on 
single particle angular momentum operators through the 
kinetic term 



i 



Define 1 = 1 — ttiluR^. In Eq. lO, if the summation can 
be replaced as an integral, then one can substitute / as 
I + mujR^ and replace the integral over / as that over 
Consequently one obtains = NmuiR^ + L'^{lo), where 
L'^{lo) is independent of lo and is thus equal to ^2(0), 
which must be vanishing. Therefore 

= NmojR^, 

which is exactly the angular momentum of a classical 
object. But is it legitimate to replace the summation 
over angular momentum eigenvalues as an integral? 

Let A6' be the range of 16*'— 6*1 in which {9'\pi\9) remains 
the same order of magnitude as {9\pi\0), while Al be the 
half width of {l\pi\l) around I = Iq. Then because of 
Eq. (jHJ, we know 

AOAi « n. 

In the absence of ODLRO, A6' < 1, thus 
Al > h. 

This allows the replacement of the summation over / as 
an integral, provided that {l\pi\l) is smooth. 

In contrast, the presence of ODLRO implies that 

A0 w 1. 

Therefore, 



Al 



which is equal to the unit difference of angular momen- 
tum eigenvalues. It is thus clear that if there is ODLRO, 
then one cannot replace the summation as an integral. 

For such a probability distribution caused by ODLRO, 
{l\pi\l) ~ Nq for / w Iq, where Nq is of the same order 
of magnitude of N, while {l\pi\l) for other values of / are 
negligible. Thus the total angular momentum is quan- 
tized as 

L, w Nolo, 

with lo determined by minimizing the Hamiltonian. 
When lli is sufficiently small, lo = 0, i.e. the system ex- 
hibits Hess-Fairbank effect. When uj is finite, lo is finite, 
but A^^o^o is less than NmuuR^. 



V. RE-DERIVATION OF SUPERFLUID 
DENSITY IN PATH INTEGRAL FORMALISM 

The analyses on solid ''He in Refs.^'^^ were based on 
an elegant path integral formulation of superfluid density 
in a rotating annulus^^ with the same geometry as in 
our consideration above. It was derived by neglecting 
the centrifugal potential. We believe that the centrifugal 
potential cannot be neglected. As this formulation of 
superfluid density is very important and widely used, it 
may be worthwhile to rederive it without neglecting the 
centrifugal potential. It turns out that it nicely remains 
the same, although the centrifugal potential is added to 
the free energy. But it seems that this is known only 
after it is checked, so it is reported here. 

We re- write the Hamiltonian in the rotating frame, al- 
ready given in Eq. as 



H 



E 



(Pj — TOV)^ 1 



2m 



- ^Nmv^ + U + V, (10) 



where, to follow RefiS^, the rotational velocity ujR is de- 
noted as V. The external potential and the interaction 
terms are schematically denoted as U and V respectively. 
U is absent in Ref.-^^, but its addition does not change 
the equations concerned. This Hamiltonian determines 
the density matrix pdm and the statistical distribution. 
One obtains^ 



P Tr{pdm) 



(11) 



where pat is the normal fluid density, P = Pj is the 
total momentum. This identity is obtained by consider- 
ing the momentum in the laboratory frame, as v is the 
container velocity in the laboratory frame. Again, note 
that the canonical momentum is the same in the labo- 
ratory and in the co-rotating frames, while it reduces to 
the kinematic momentum in the laboratory frame. 
Because 



dH 

dv ' 



Eq. Hll|l can be re-written as 



PN dF 
— Nmv = - — , 
p ov 



(12) 



where F = —kTlii[Tr{pdm)] is the free energy in the 
co-rotating frame. Therefore the superfluid fraction is 



PS 
P 



1 



d{hniv'^ 



Thus the free-energy change due to the rotation of the 
container, up to the order of w^, is 



mv^ PS 
2 ^ p 



1), 



(13) 
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from which it can be confirmed that the centrifugal po- 
tential mv^ /2 indeed cannot be ignored, since it is no less 
than the other term {ps / p)mv'^ /2. 
In the path integral calculation, 

g-/3A_F _ /Pdm(X,X;/3;v)dX 

/ prf„(X,X;/3;v = 0)dX' 

where X represents the configuration of the particles. 
The "gauge term" — mv in the kinetic energy term can 
be transformed away, by adding, in the density matrix 
elements, a phase factor in winding the periodic system, 
like in Eqs. © and {Tj). Consequently, one can replace 
p<im(X,X;/3; v) as the density matrix pd,„(X, X; /3; v) 
corresponding to the Hamiltonian without the "gauge 
term" , while multiplying it by a phase factor due to the 
total paths WL of the N particles winding around the 
system, where W is the winding numberi22i 

= exp(-/3fl^), where H = H{v = 0) - Nmv^/2. 
H{v = 0) is just in Sec. III. We obtain 

where the average that of the density matrix with v = 0. 
Consequently, up to the order of , we have 

^^^^W m(I^_ 

which, together with Eq. ((T^ . yields 

PS _ m{W^)L^ 
~p ~ 3f3h^N ' 

which is the same as that given in Ref.^2, They remain 
the same even if v is not a small quantity, for the reason 
is that V is independent of the particle configuration. 

VI. NO ODLRO IN NOSANOW-JASTROW 
WAVEFUNCTIONS 

Now we turn our attention to solid ^He. For a commen- 
surate solid at rest, each atom occupies a lattice site. Be- 
cause of quantum mechanical zero-point motion, which 
is large in solid helium, around the neighborhood of each 
lattice site, there is a finite region in which the wave- 
function is nonvanishing. With the exchange effect put 
aside first, the wavefunction is localized around each lat- 
tice site, i.e. it decays from the maximum at the lattice 
site. Let's denote the wavepacket of atom i as w{ri — Qi), 
where is the actual position of the atom, represents 
a lattice site fixed in the solid. The Hartree approxima- 
tion of the wavefunction of the solid helium is the product 
of these single-atom wavefunctions, i.e., 

JV 

<^H = l[w{r,-Q,), (14) 

i=l 



which was indeed used in the earliest (unsatisfactory) 
variational calculations of solid heliumiS Later works, 
starting by NosanowfS& took into account the two-particle 
short-range correlation by multiplying the Hartree wave- 
function by the Jastrow factor. 

To account for the exchange effect due to overlap be- 
tween neighboring single-particle wavepackets, one also 
needs to consider the wavefunction symmetrized over all 
the atoms; the detailed nature of the exchange effect is 
then determined by the Hamiltonian. With symmetriza- 
tion, the Hartree approximation is improved to Hartree- 
Fock approximation, 

1 " 

p i=i 

where P represents TV! permutations of the N lattice 
sites {Qi}. The symmetrization can be made on either 
the particle positions {r^} or the lattice sites {Qi}. We 
choose the latter for easier manipulation below. 

The symmetrized Nosanow- Jastrow wavefunction is 

N 

'^SNJ = KY,PllHr,-Q,)l[l[fjk, (16) 

P i=l k j<k 

where K is the normalization constant, 

fjk = /(-u(|rj -rfcl)) 

is the Jastrow (or, to be historically precise, Bijl-Dingle- 
Jastrow) function. /(— u(r)) attains a maximum larger 
than 1 at a certain distance rg , and it is constrained to be 
/ ^ as r — > 0, and / ^lasr— >ooorr>cr where a 
is a parameter. Note that Yik Y[j<k fjk is automatically 
symmetric for all particles. 

Our consideration is about a thin cylindrical bulk, 
Ti — {R, R9i, Zi). Especially, the periodic boundary con- 
dition in coordinate 9 should be taken into account in an 
essential way. It implies that the Wannier-like function 
w must be of the formSi 

oo 

w{r-Q) = A J2 u'(r-Q-7G), (17) 

7— — oo 

where A is the normalization constant, G ~ 2ttR9 rep- 
resents the circumference, 7 represents integers, w is the 
(real) Wannier-like function for the infinite interval. Each 
w extends over a finite range, much smaller than the sys- 
tem size, but finite overlap is allowed. ±00 in the sum- 
mation can be understood as two bounds which can be 
arbitrarily large. Thus 

w(r)iD(r - S) w w^{r) exp(-|S|/c), (18) 

where S is an arbitrary vector, c is a length scale less 
than the lattice constant. Consequently, the normaliza- 
tion constant A in Eq. ((T7jl is A « (X]7 7' 6xp(— 17 — 
7'|G/c))-i/2. 
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Moreover, it can be found that 

w(r)w(r-S) < w^(r)exp(-|5e|/c), (19) 

where Sg is the 9 component of S modulo ±G such that 
|<5e| < G/2, i.e., \Sg\ is the shortest ^-component of the 
distance between the two physical points represented by 
r and r — S . 

We now set out to show that there is no ODLRO or 
BEC in or ^hf or ^snj, by examining the one- 
particle reduced density matrix 

pi(r,r') = N j dY2- • • drA,$(r, rz, • • • , rAr)$(r', ra, • • • ,rAr) 

for the ground state wavefunction $ of the form of 
or <^HF or $SArj. 

Though trivial, it is instructive to first consider . It 
is straightforward to integrate out r2, • • • , tat, and obtain 
/rji(r,r') = Nw{y - Qi)w(r' - Qi), for which Eq. ^ 
directly leads to 

/9i(r,r') < iViD2(r)exp(-|a;-2:'|/c), (20) 

where x = R6 denotes the 6'-component of r. Of course, 
w^(r) < 1. Thus yOi(r, r') — > as |a; — x'l approaches the 
system size, i.e. there is no ODLRO or BEC in . 

Now consider the Hartree-Fock wavefunction p . In 
the expansion of pi, suppose the lattice sites in the first 
$ are denoted as {Qi} while those in the second $ are 
denoted as {Q-}. The exponential decay of the overlap 



where /{j = f{—u{\r[ — r^D). The fraction factor in 123(1 
must be bounded by a finite number. Clearly, pi(r, r') 
tends to exponentially vanish as |x — approaches the 
system size. Thus there is no ODLRO or BEC in ^snj 
either. It can be seen that our argument is not disrupted 
by the thermodynamic limit N oo. 

Furthermore, the argument can be straightforwardly 
generalized to a finite temperature, in which each energy 
eigenfunction is of the form of or ^hf or ^snj- The 
finite-temperature density matrix is the thermal average 
of the density matrices of the eigenfunctions. For an 
infinite sample, w would simply be w, the conclusion of 
no ODLRO can still be obtained, in a simpler way. 

Therefore, although there is ODLRO or BEC in the 
Jastrow wavefunction alone, which describe liquid he- 
lium,^-!, they are dominated by the localized one-particle 



between single- particle wavefunctions, Eq. H19() . implies 
that among the (A^!)^ terms in the expansion of pi, one 
can neglect each term in which Qi 7^ Q^ for at least one 
of « = 2, • • • , A^. Consequently, there arc only A^! remain- 
ing terms, in each of which Q^ — Q'^ for i — 1, - ■ ■ ,N, 
then r2, • • • , rjv are subsequently all integrated out. This 
iV! is cancelled by the A^! in the normalization constant. 
Hence, for large jx — pi(r,r') for ^hf behaves in 
the same way as for the Hartree wavefunction, given in 
Eq. 1(201). This proves there is no ODLRO or BEC in 
^HF either. 

The argument can be extended to symmetrized 
Nosanow- Jastrow wavefunction ^snJi which can be re- 
written as 

N 

'^SNJ^KJ2Pl[Hr^-Q^)Y[fJ^], (21) 
P i=l j<i 

where P represents the permutation of the A' lattice 
sites {Qi}. Ylj^i fji is a function of ri, • • • , ri, and re- 
duces to 1 for i — I. For each term in the expan- 
sion of pi, consider w(ri - QOw(r.j - Q'i)illj<i fjt)'^ < 
w'^{r)exp{-\Q^g - Q'^g\/ €)([{. ^^fji)"^, where Q,e is the 
9 component of Qi. It can be seen that the short-range 
Jastrow factor does not change the nature of long-range 
exponential decay. Therefore, the cross terms, in which 
Qi 7^ Q^ for at least one of i = 1, • • • , A^, exponentially 
decay, and are negligible in comparison with the remain- 
ing terms. Consequently, 



(22) 
(23) 



wavefunctions. This is a difference between liquid and 
solid. The argument extends that of Penrose an Onsager 
about no BEC in a solid^ to the case with zero-point mo- 
tion, exchange effect, as well as short-range correlation. 



VII. NO SUPERSOLIDITY IN 
NOSANOW-JASTROW WAVEFUNCTIONS 

As ODLRO is a sufficient condition of NCRI, it is not 
redundant to demonstrate that there is no NCRI in 
or ^HF or ^SNJ, as we now explicitly do in the follow- 
ing. We adapt the method of Kohn used in discussing 
electronic insulating statei^ 

Recall that the eigenfunctions and energy spectrum is 
determined by Hq, as in Eq. (|SJ|. The idea is the foUow- 



pi(r,r') 



Nw{r^ - Q^)w{r[ - Qi) n.>i / ^'^(i-. - QO(ni<,<. 
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ing. For every eigenfunction \E'q(w = 0) of Ho{uj = 0), 
where a is the index for different eigenfunctions, be it of 
the form of or ^hf or ^snj, we show that there is a 
corresponding eigenfunction ^'^(w ^ 0) of Ho{uj ^ 0), 
and that its eigenvalue remains the same as that of 
Ho{u> = 0) for = 0). 

Ho{llj 7^ 0) is simply related to Ho{uj = 0) by a "gauge" 
transformation, but one should be cautioned by the re- 
quirement of the SVBCf^iSi Eq. In an infinite inter- 
val, for a localized single-particle eigenfunction w{r) of 
a single-particle Hamiltonian, the correct eigenfunction 
wavefunction for lu ^ is 

_ _ imujx 

w [uj; r) ~ w[r) exp( — - — ), 

where, as above, x — R9. 

Therefore, for a many-particle eigenfunction ^a{^ = 
0) of Hoiu) = 0), given by or ^hf or ^snJi one may 
construct the corresponding eigenfunction <&ct(i^ 7^ 0) of 
Hq[lo 7^ 0) in a similar way, by replacing every single- 
particle ^(r) as ui'(r). The presence of Jastrow factor 
does not affect this. 

On the other hand, by using Eq. '^a{'^ = 0) can 
be written as 

oo 

*„(c. = 0) = A^ ^ ^r({ra), 

r=-oo 

where is obtained from $ by shifting the centers of 
the single-particle wavepackets w from {Qi} to {Qi + 
7iG}, with X]i7j — here !> is of the form of = 
a=i ^(r. - QO, or ^HF = ^ Ep ^nti " QO, 

or ^SNJ = Ep ^nti ^(r. - Q.) Yik Y\j<k /.fc- 

Following the argument in Refi^, using the exponen- 
tial decay of the overlap as given in Eq. (fH^ . and very 
similar to the argument in last section, it can be shown 
that for r ^ r' and arbitrary a and a', <i>Q,r and ^a',T' 
have exponentially vanishing overlap and give vanishing 
matrix element for an arbitrary one-particle position op- 
erator. 

Consequently, it can be found that the corresponding 
eigenfunction of i?o(w), satisfying the SVBC, is 

vI/„M= £ <i>„,r({rj)exp[^(^^^,-27rr)]. 

r=-oo j 

(24) 

Because of exponentially vanishing overlap between '^a,v 
with different values of F, it is clear that 

with 

E^{lo) = E^iio = 0). 

It is thus proved that every eigenvalue Ea{Lu) of Hq{uj) 
is independent ofuj. Interestingly, the argument has gone 
through even in presence of the Jastrow factors. 



In fact, the explicit construction of the wavefunction 
here confirms the principle, established by Leggettj^ that 
the system is non-superfluid if for the wavefunction of the 
rotating system, the SVBC can still be satisfied with- 
out causing the energy to be increased by the rotation. 
Leggett already applied this principle to the Hartree-Fock 
wavefunction. 

Therefore, for a commensurate quantum solid de- 
scribed by Hartree or Hartree-Fock or Nosanow- Jastrow 
wavefunction, even though the exchange effect, large 
zero-point motion and short-range correlation are taken 
into account, the free energy is of the form of Eq. 
with Fq independent of uj. This indicates that it cannot 
a supersohd. 

In our argument, the localized single-particle wave- 
functions play a crucial role. Obviously, the situation 
would be different when there exist vacancies or intersti- 
tials or both, which makes the wavefunctions extended. 
The recent experimental result of Kim and Chan poses a 
significant challenge. The difficulty might be resolved if 
an extended factor is found in the actual wavefunction. 



VIII. SUMMARY 

To summarize, we have offered some analytic argu- 
ments concerning the existence or non-existence of su- 
perfluidity or supersolidity behavior. This work might 
be useful for further investigations on the cause of su- 
persolidity. It might be helpful in supplementing the 
understanding of the relevant classic literature, and in 
clarifying which and which specific features are counter- 
parts between superfluidity and superconductivity. 

Our argument seems to suggest that ODLRO is indeed 
generically sufficient for superfluidity even in a finite tem- 
perature, a question which seems to have remained not 
entirely resolved previously. 

Our discussions start with a synopsis, in Sec. II, on the 
Hamiltonians and the free energies in the co-rotating and 
the laboratory reference frames, as well as their relations 
with rotational inertial and superfluidity density. 

In Sees. Ill and IV, from the presence of ODLRO, we 
make a general derivation of the most basic manifesta- 
tion of superfluidity, namely the Hess-Fairbank effect or 
NCRI, i.e., the refusal of the Bose system to follow the 
rotation of the container, by using a method of Yang in 
treating superconductivity in a magnetic field. We also 
derive the quantization of angular momentum as a con- 
sequence of ODLRO, by borrowing a method of Bloch 
in studying superconducting persistent current. In Sec. 
V, we rederive the path integral formulation of the super- 
fluid density without neglecting the centrifugal potential. 

In Sees. VI and VII, we consider the variational wave- 
functions which have been used in solid helium calcu- 
lations, namely, the Hartree, the Hartree-Fock and the 
symmetrized Nosanow- Jastrow wavefunctions. The non- 
superfluidity in the Hartree-Fock wavefunctions was al- 
ready noted by Leggett from its disconnectivity^S. We 
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show that there is no ODLRO in these trial wavefunc- 
tions, for both an infinite sample and that confined in 
a cylindrical annulus. Moreover, by extending a method 
originally due to Kohn in discussing electronic insulating 
states, we explicitly demonstrate that there is no NCRI 
behavior in a commensurate quantum solid described by 
those trial wavefunctions, even if there exist large zero- 
point motion, finite overlap between wavepackets and ex- 
change effect. In this argument, the constraint of SVBC 
in the angular direction is carefully taken into account. 
The explicit construction of the wavefunction under the 
rotation is consistent with the early arguments of Leggett 



in terms of the connectivity properties^ii^ 
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Erratum 

[published as Phys. Rev. B 74, 029901 (E) (2006)] 



In Sees. VI and VII, the treatment of symmetrization 
was inadequate. The problem is fixed below. Our condi- 
tion is more relaxed: the conclusion holds as far as the 
overlap between two atomic Wannier functions decays 
exponentially, with the decay width c much less than 
the system size, rather than the lattice constant. That 
is, M)(r)w(r — S) « iD-^(r) exp[— |S|/c], for the Wannier 
function w in an infinite interval, or 'w(r)'w{r — S) < 



'it;^(r) exp[— |0(S)|/c] for the Wannier function w in the 
annular geometry, where 8(S) denotes the shortest 9- 
component of S. As in the original paper, we discuss the 
annular geometry; the case of infinite interval is simpler. 

Consider the Hatree-Fock wavefunction ^hf — 
Z3 Hi^i — QPi), where P represents N\ per- 
mutations of the N lattice sites {Qi}, {Pi} are the in- 
dices for the N lattice sites in permutation P. It can 
be found that D = 1/N\VA, A = J^pUiOiP,, where 
OiPi = J driw{ri — Qi)w{ri — QpJ. The one-particle 
density matrix is thus found to be 



Pi (r, i") = ^ w;(r - Qpju;(r' - Qp; ) H / dr^w{r, - Qp^)w{ri - Qp- 

p.P' 1^1 



^ WxT. Mr-Qp,)Mr' - Qp;)nexph|e(Qp. - Qp;)|/c] 



iV!A 



p.P' 



^ Ma ^ ~ QpJ cxp(-|e(r - r' - Qp, + Qp;)|/c) exph|e(Qp, - Qp.)|/c] 
p,p' 
N -Nl 



< 



■exp[-|e(r-r')|/c]. 



In the derivation, the relation ^ ■ Qp. — J2i Qf/ has 
been used. Therefore, pi{r,r') ^ as |8(r — r')| ap- 
proaches the system size. The argument is valid in the 
thermodynamic limit, as A is of the same order of mag- 
nitude of A^!. 

Now consider the symmetrized Nosanow-Jastrow 

wavefunction ^snj = KNY.pI\^=i'^i''^t ~ 
QPi)Ilj<kfik, where fjk = f{\rj - rfe|) is the Jas- 
trow function, K]^'^ = J2p,p' Hi I ^ri Jlj - 
Qp,)w(r^ - Qp/)(nj<fe /jfe)^- In a way simi- 
lar to the above derivation for ^hFi it is found 
that pi(r,r') < NK%Y.P,P'Wi^ - QPi)w{r' - 
Qp;)exp[-|e(Qp, - QpO\/c]Ui^Jdvilli^,f{\r - 
r,|)/(|r' - r,|K(r, - Qp^mi<j<k fjk)' < 
Arexp[-|e(r - r')\/c]F^K%j:pp,U.^, J dnw^ir, - 
Qfi)(ni<7<fc /jfe)^' where F is the maximum value of 



the Jastrow function. It can be seen that pi(r,r') 
as |0(r — r')| approaches the system size, and that the 
argument is valid in the thermodynamic limit. 

Thus it is proved that there is no ODLRO in ^hf 
or $swj under the condition stated above. The sym- 
metrization should be dealt with in a similar way in 
Sec. VII, where it is shown that ^a,r and ^a'.r', with 
T ^ T', have exponentially vanishing overlap and give 
vanishing matrix element for an arbitrary one-particle 
position operator, which leads to the conclusion of no 
supersolidity in these wavefunctions. 

As noted in Ref. 1, it had been shown previously that 
there is no ODLRO in Hatree-Fock wavefunctions under 
the condition that the sum of the overlap integrals be- 
tween any site and its neighbors is less than unity^i. 

I thank Xin-Cheng Xie for questioning the original 
treatment of the symmetrization. 
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